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Laser Interferometer Space Antenna (LISA) will routinely observe coalescences of supermassive 
black hole (BH) binaries up to very high redshifts. LISA can measure mass parameters of such 
coalescences to a relative accuracy of 10"*-10"^ for sources at a distance of 3 Gpc. The problem of 
parameter estimation of massive nonspinning binary black holes using post-Newtonian (PN) phasing 
formula is studied in the context of LISA. Specifically, the performance of the 3.5PN templates is 
contrasted against its 2PN counterpart using a waveform which is averaged over the LISA pattern 
functions. The improvement due to the higher order corrections to the phasing formula is examined 
by calculating the errors in the estimation of mass parameters at each order. The estimation of the 
mass parameters M and rf are significantly enhanced by using the 3.5PN waveform instead of the 
\^ , 2PN one. For an equal mass binary of 2 x IQ^Mq at a luminosity distance of 3 Gpc, the improvement 

• in chirp mass is ~ 11% and that of is ~ 39%. Estimation of coalescence time tc worsens by 43%. 

' The improvement is larger for the unequal mass binary mergers. These results are compared to 

^N) ^ the ones obtained using a non-pattern averaged waveform. The errors depend very much on the 

location and orientation of the source and general conclusions cannot be drawn without performing 
Monte Carlo simulations. Finally the effect of the choice of the lower frequency cut-off for LISA on 
the parameter estimation is studied. 
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I. INTRODUCTION 



(N 
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\ Most of our near-by galaxies harbour supermassive black holes (BH) at their centre If this is the case, then 
merger of such galaxies would produce a binary system composed of two supermassive BHs. Simulations indicate that 
there could be primordial supermassive BH binaries at the centres of the first galaxies An understanding of the 
formation and evolution of these binaries is very important from the view point of cosmology and structure formation. 
Many of these binaries coalesce under gravitational wave (GW) radiation reaction within Hubble time. Supermassive 
'f , BH binaries in the mass range 1O^-1O^M0 will emit gravitational waves (GW) of frequency 10~^-10~^ Hz during its 
J~^_ ■ adiabatic inspiral phase which can be observed by the proposed space-borne GW missions such as LISA H with 
rrl high ( ~ a few thousands) signal to noise ratio (SNR) up to very high redshifts (~ 10). 

K>" ■ Different authors have investigated the implications of these observations in the context of astrophysics, cosmology 
k> I and testing general relativity and its alternatives. LISA observations of BH coalescences can be used to study the 
'jjj . growth of BHs as the universe evolved and for mapping the distribution of BHs as function of the redshift 0, |E 01 ■ 
' LISA will be able to measure luminosity distances to the sources with an accuracy ~ 1 — 10%. If the redshift associated 
with the event is known by electromagnetic observations, these sources can be used as very high precision standard 
candles and to study the distance- redshift relation 01 ■ Ref ^3 discussed the potential of LISA to observe binaries 
containing a BH in the intermediate mass regime (~ IO^Mq) and use it as a probe of strong field aspects of gravity. 

LISA could probe many strong gravitational field effects which are not possible to explore with other observational 
means. Both inspiral and ring-down GW signals can be used for this. Refs 0,^3 studied the possibility of using the 
quasi-normal mode oscillations to test the no-hair theorem of general relativity since these modes will be characterized 
only by the mass and angular momentum of the BH (in general relativity). Hughes and Menou examined another 
possibility ^3 if LISA detects both inspiral and ring-down signals from the same source. By measuring the total 
mass independently from both the signals one can estimate the mass difference which will be the mass-energy lost 
due to GWs. They suggested that an extension of this idea including spin effects could in principle test the BH area 
theorem. Further, inspiral of a stellar mass BH into a SMBH will be another interesting source for LISA using which 
many properties of the central SMBH can be probed including the possibility to map the spacetime following the 
geodesies of the stellar mass BH (See e.g.pj, lli) and measuring the multipolc moments of the spacetime. 

LISA will provide an unique opportunity to test general relativity and its alternatives. Will and his co-workers have 
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discussed the potential of LISA to test general relativity as well as its alternatives like Brans-Dicke theory and massive 
graviton theorie s |l6l Il7l | . Recently this issue was discussed in a more realistic scenario of spinning binaries using 
2PN phasing 0, Q . Blanchet and Sathyaprakash proposed another test based on post-Newtonian (PN) GW phasing 
formula by measuring the 1.5PN GW tail effect and showing how it can be used as a test of general relativity [ia.[l9l|. 
This proposal was recently generalized to higher order terms in the phasing formula by Arun et al pol | and it was 
argued that such a test would allow one to probe the nonlinear structure of gravity |2ll |. 

A very accurate parameter extraction scheme is central to performing all these analyses. A parameter estimation 
scheme based on matched filtering, similar to that for the ground based detectors such as LIGO and VIRGO, will be 
employed for LISA also. An efficient matched filtering would in turn demand a very accurate model of the gravitational 
waveform. In order to compute the gravitational waveform from a compact binary system, one solves the two-body 
problem in general relativity perturbatively using different approximation schemes since no exact solutions for this 
problem exist till date. The final waveform can be expressed as a post-Newtonian expansion which is a power series 
in v/c where v is the gauge independent velocity parameter characterising the source (See Ref. [2^ for an exhaustive 
review on the formalism). In our notation - refers to half a PN order. 

Since the information about the phase is more important for the process of matched filtering, one uses a simplified 
model of the inspiral waveform (the so called restricted waveform) where phase is modelled to a high PN order retaining 
the Newtonian amplitude. In doing so, one is neglecting the effect of other harmonics |2^ l24j in the amplitude of 
the wave and also the higher order PN corrections to the dominant harmonic at twice the orbital frequency. In 
the present study, we deal only with the restricted waveform in the Fourier domain obtained using stationary phase 
approximation. The phasing formula for nonspinning binaries, is presently complete up to 3.5PN order p5ll26ll2^ l28l | . 

The implications of the higher PN order phasing in the context of parameter estimation problem has been investi- 
gated by different authors. Based on the framework set up by Refs [23.l30|. Cutler and Flanagan [sj investigated the 
importance of the 1.5PN phasing formula [s^- The effect of including spin-orbit coupling parameter at 1.5PN order 
into the space of parameters was one of interesting issues addressed. Two independent works by Krolak et 
and Poisson and Will 01 analysed the problem of parameter estimation using the 2PN phasing formula of Ref. 
Inclusion of the spin-spin coupling term at 2PN and its effect on the errors of other parameters was the focus of their 
analysis. 

Recently Arun et al. investigated the effect of the 2.5, 3 and 3.5PN terms for the parameter estimation of 
nonspinning binaries (a similar work was carried out independently by Berti and Buonanno |3^). Using covariance 
matrix calculations, they inferred that by employing the 3.5PN phasing instead of the 2PN one, the improvement 
in estimation of errors in chirp mass and symmetric mass ratio can be as high as 19% and 52% respectively for the 
ground based detectors such as LIGO and VIRGO. 

Cutler was one of the first to address the problem of parameter estimation in the LISA context [s^- He used 
the 1.5PN waveform including the spin-orbit effect and studied the estimation of errors associated with the mass 
parameters as well as distance and angular resolution of the binary. Seto investigated the effect of finite arm length 
of LISA using 1.5PN phasing [s^- Vecchio revisited the parameter estimation with the 1.5PN waveform [s^ where 
he used the waveform for circular orbit but with "simple precession" (as opposed to the non-precessing case of [s^ ) 
and examined the implications of it for the estimation of distance and angular resolution. Various aspects of the 2PN 
parameter estimation, such as the spin-spin coupling, was investigated by different authors ISJ3L -R-^f studied 
the effect of spin terms in testing alternate theories of gravity with the LISA observations. Refs |a, Q also addressed 
the issue of mapping the merger history of massive BHs using LISA observations in the 2PN context. While all 
these calculations are within the restricted waveform approximation where the PN corrections to the amplitude is 
completely neglected, there have been investigations about the effect of including these amplitude corrections in the 
context of parameter estimation EdI • 

Other than the covariance matrix approach, which is valid only in the high SNR limit, there have been proposals 
in literature addressing the parameter estimation problem using Monte-Carlo methods. In Ref the authors 
compared the error estimates obtained using the covariance matrix with the Monte-Carlo simulations. Recently 
parameter estimation schemes based on Bayesian statistics using Markov chain Monte-Carlo (MCMC) method also 
has been proposed and implemented in the ground based detector context ji^ and the LISA case [47l l48l |. 
Using 2PN waveform Ref |47l | found that posterior parameter estimation distribution of the extrinsic parameters 
obtained using MCMC methods are in excellent agreement with those computed using Fisher matrix whereas there 
is a systematic overestimate by Fisher matrix for the intrinsic parameters. 

In the present work we extend using the covariance matrix formalism, the parameter estimation problem for 
supcrmassive BH binary inspirals in LISA and study the implications of higher PN order terms. Coalescences of 
BHs of masses 10* — 10^ Mq at a luminosity distance of 3 Gpc are considered. We assume LISA will observe these 
events for one year duration. Using the 3.5PN phasing we calculate the errors associated with the estimation of the 
mass parameters and coalescence time and compare them with the corresponding 2PN results. We also study, for 
chosen source locations and orientations, the effect of orbital motion for parameter estimation by comparing these 
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results with the other two cases: one where LISA is assumed to be a single interferometer and another when LISA is 
considered to be a two detector network. Conclusions drawn from the exercise above is far from general and have to 
be supplemented in future by rigorous Monte Carlo simulations. 

The rest of the paper is organized as follows. Sec. ITTl discusses all the necessary inputs required for the paper such as 
a brief introduction to parameter estimation using covariance matrix, noise model for LISA, model for the waveform 
and some other conventions followed in the paper. Sees IIIII and IIVI discuss the main results and their implications 
and Sec. Ivl provides the summary and future directions. 



II. PARAMETER ESTIMATION FOR LISA WITH 3.5PN PHASING 



A. Parameter estimation using the covariance matrix 



We summarize the theory of parameter estimation in the context of Gaussian random detector noise, addressed 
in the GW context first by Finn and Chernoff js^l and implemented in detail by Cutler and Flanagan [sj ■ Let 
us assume an inspiral GW signal is detected meeting the necessary detection criteria and one needs to extract the 
intrinsic and extrinsic parameters from the signal by matched filtering. 

For sources like inspiralling compact binaries, where a prior source modelling is possible to predict the gravitational 
waveform, matched filtering is an ideal method both for detection as well as parameter estimation of the signal |49l | . In 
matched filtering, the detector output is filtered using a bank of theoretical templates with different signal parameters. 
The parameters of the template which obtains the best signal-to-noise ratio (SNR) gives the "measured" values of the 
signal parameters. These values, in general, will be different from the "actual" values due to the presence of noise. 
The problem of parameter estimation addresses the question of how close are the measured values to the actual ones 
and what the associated errors are in the estimation of different parameters. For a given signal, different realizations 
of noise lead to different sets of best-fit parameters of the signal. When the background noise is a stationary, random, 
Gaussian process, at high enough SNRs the best-fit values of the parameters have a Gaussian distribution centered 
around the actual values of the parameters. 

If A* denotes the actual value of the parameters and A'+AA', the measured value, then the root mean square 
difference AA* obeys a Gaussian distribution: p(AA*) cx exp (— F,,, AA' AA-'/ 2) where F^, the Fisher information 
matrix constructed from the Fourier domain representation of the waveform, is given by 

Here, hi{f) :— dh{f)/dy, h{f) is the Fourier domain gravitational waveform and Sh{f) is the (one-sided) noise power 
spectral density of the detector, /in and /gn denote the lower and upper limits of integration which is defined later 
in the section. It also follows that the root-mean-square errors arc given by cr' — a/S", where E = F~^ is called 
the covariance matrix. The non-diagonal elements of the covariance matrix are related to the correlation coefficient 
between two parameters A* and A^ : c*^ :~ Repeated indices are not summed over in the above expressions. 

Finally the SNR can be expressed in terms of the Fourier domain signal h{f) as 

In the present case the A' denoting our chosen set of parameters are given by {tc, </>c, A^, ??, D]^, (psi (f'L, (^Si ^'l}- In ths 
case where pattern functions are not included the above set reduces to {tc, 0c,.A/^, The additional elements of the 
parameter set denotes the distance, orientations and locations of the source in the sky specified with respect to the 
fixed solar system based coordinate system. 

In the above integrals, the upper limit of integration is /fin = Min[/iso, /end], where /iso is the frequency of the 
innermost stable circular orbit for the test particle case, fiso — (6"^/^7rm)~^ and /upper corresponds to the upper 
cut-off of the LISA noise curve /end = IHz. We have chosen the lower limit of frequency /in = Max[/in, /lower] where 
/in is calculated by assuming the signal to last for one year in the LISA sensitivity band and /lower, the low frequency 
cut-off for LISA noise curve, is assumed to be lO^^Hz^. This is a rather optimistic choice for the lower cut-off for 



^ Another way of choosing the limits of integration is to calculate the time over which the signal will last once it enters the LISA band. 
See 1^ for example, where the duration of the signal is calculated using the expression for It assumes /lowcr — 

lO-'Hz and 3 year 

mission time for LISA. 
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LISA and in Sec. IIII d wc estimate its effect on parameter estimation by comparing it to the more modest choice of 
10--^ Hz. 

We follow the noise model of LISA as given in Sec. IIC of Ref which is a slightly modified version of The 
noise spectral density consists of a non-sky averaged part [sof and confusion noise due to the galactic and extra galactic 
white dwarf binaries |51i, iS^j . The explicit expressions for these components can be found in Eqs (2.28)-(2.32) of 
Ref 0. 



B. LISA detector configuration and the waveform model 

LISA is a three-arm interferometer where each arm has a length of 5 x 10^ km. These arms form an equilateral 
triangle and move in a heliocentric orbit with a 20° lag to the earth and the plane of the detector tilted at 60° with 
respect to the ecliptic (See 0, 113 for details). LISA with three arms is essentially equivalent to a pair of two-arm 
detectors capable of simultaneously measuring the two polarizations of the incoming gravitational wave. We consider 
two cases: one where we assume the estimation of mass parameters is not affected because of their correlations with the 
angular variables and second when we estimate the associated errors with angular variables and luminosity distance 
of the source. Since the information about the angular variables are encoded in the so-called pattern functions which 
describes the orbital motion of LISA, in the first case we use a waveform which is averaged over the pattern functions. 
In the second case, we do not average over the pattern functions and use the information from the LISA orbital motion 
to discuss the estimation of angular resolution and luminosity distance to the source. Further, in the second case we 
consider cases when (i) LISA is a single two arm-detector and (ii) as a two detector network in order to understand the 
effect of network configuration for parameter estimation. The LISA antenna patterns, describing its orbital motion, 
is given in js^l which is used for the present study (also see Appendix A of Ref 0| for these expressions). 

Unlike the ground based detectors where the two arms have an angle 90°, the LISA arms are at 60°. As shown 
by Cutler [s^], the relative strain amplitude (which is the gravitational waveform) in the LISA case can be simply 
be related to the 90° interferometer case by multiplying the latter by a factor •\/3/2. Using this input, the Fourier 
domain waveform within the stationary phase approximation can be written down as 0, Is^ 

hM) = ^Ar"^e'^'^'\ a = 1,11, (2.3a) 

where a labels the interferometer, / the GW frequency and M., the chirp mass which is related to total mass 
m = TOi -|- 777.2 and symmetric mass ratio t] = 77117772/777^ hy M = rf^^ m. The luminosit y d istance to the source is 
denoted by D^. The GW phase 7/) (f) appearing in the formula is completed up to 3.5PN |27ll28l l53 | and its Fourier 
domain representation is given in |5J, |55j . We find it more convenient to write it as 

fe=7 

i:U)^2^ft,~^, + Y,akv\ (2.4) 

fe=0 

where v = {nm fY^^ is the PN variable which is related to gauge independent source velocity in system of units where 
G = 1 = c which we follow henceforth in the paper. Eq (3.4) of gives the ak for different values of fc = • • • 7. 
In the case where we do not average the pattern functions, the waveform can be written as 

where ipp^a{t{f)) and LpD{t{f) are the polarization phase and Doppler phase respectively [s^ . Aa{t{J)) correspond 
to the amplitude modulations induced by the EISA's orbital motion. Aa{t{f)) and fp^aitif)) depends on the pattern 
functions F"{t) and F"{t) and hence vary with time. For the explicit expressions for ipp,a and ipjj, we refer the 
readers to Rcfs 0,113. For 3.5PN accurate expression for t{f) we use the following relation 

2^tif) = '-fl. (2.6) 

This can be rewritten as 

7 

t{f) = (2.7) 

k=0 
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and values of is given in Refs |5J,l55j which can readily be used. 

For calculations where LISA is assumed to be a two detector network, we calculate the SNR and Fisher matrix 
using 



P""™- Vp?+P?i, (2.8) 

T^Nctwork T^I I T^II /o n\ 

i q6 —'^ab+^ab- l^-yj 

The errors for the two detector case are obtained inverting the total Fisher matrix following the procedure outlined 
in Sec. ESI 

Throughout the paper we assume a cosmological model with zero spatial curvature fi^ = 0, CI a + = 1 and 
Hubblc's constant to be Hq = 70 km s~^ Mpc~^. fl\ and CIm refers to the contributions to the total density from 
matter and cosmological constant. The luminosity distance is given by 

D.^iiLff 'ISl (MO, 

Ho Jo [fiMil + z'^ + nAf^ 

where z denotes the redshift of the source. 

We calculate the Fisher matrix for the different configurations of LISA using the corresponding waveforms and 
invert it to get the covariance matrix. The elements of the covariance matrix are used for discussing the errors and 
correlation coefficients of different parameters in the next section. While discussing the trends with the FN orders, 
its useful to keep in mind that when pattern functions are included, there are additional FN expansions coming from 
t{f) both in the amplitude and phase, apart from the usual phase ipiDi which can influence the results. 



III. PARAMETER ESTIMATION WITH PATTERN AVERAGED WAVEFORM 

In this section we discuss the parameter estimation in the LISA case using pattern averaged waveform of Eq (|2.3|l . 
Parameter estimation with the non-pattern averaged waveform (as in Eg I2.5|l is more complex since the estimated 
errors strongly depend on the location and orientation of the source (see the Sec. lIVI below) which enter the calculation 
via the pattern functions. The ideal way to deal with the situation will be to perform Monte Carlo simulations for 
different binaries located and oriented randomly in the sky 0, S S US ■ A recent Monte Carlo simulation (6| , which 
addressed the parameter estimation problem using the 2PN waveform including spin effects, compared the results of 
their simulation with the result they obtained using pattern averaged waveform (see Tables V, VI and VIII of Ref j^). 
They found that the results in both the cases are in excellent agreement. The results presented in this section about 
the improved parameter estimation with the pattern averaged waveform, may hence give a reasonably good idea about 
the full problem where the LISA pattern functions are included and it is considered to be a two detector network. 
We emphasize that the results quoted here have to be supported by Monte Carlo simulations similar to Ref (see 
concluding remarks in Sec. Iv|) . 



A. Improved parameter estimation of equal mass binaries with the 3.5PN phasing 

We discuss the performance of the 3.5PN restricted waveform from the parameter estimation point of view for the 
pattern-averaged case discussed above. Our aim is to study the variation of errors in different parameters with the 
total observed mass^ of the binary for different FN orders. This would not only give us an idea of the improvement 
brought in by the use of higher order phasing but also about the convergence of the FN series for the problem of 
parameter estimation. We have checked our codes by reproducing the results at 2FN with that in Table III and V 
of Ref for the nonspinning case. The important results of our study are discussed in detail in what follows. The 
errors in estimation of different parameters for a 2 x 10® A/© binary at 3 Gpc is provided in Table I for different FN 
orders in the phasing. 



^ By total mass, we always refer to the total rcdsliifted mass m'(l -|- z), where m' is the actual source mass and z is the redshift of the 
source. This is the mass parameter that is measured by the GW observations. 
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FIG. 1: Variation of errors with the total observed mass for different PN order restricted waveforms for LISA. Pattern averaged 
waveform is used. The convergence of the results is evident in both the cases. Sources are assumed to be at 3 Gpc. 



B. Improvement in estimation of mass parameters and PN convergence 

Improvement due to higher order terms: We plot in Fig. ^ the variation with mass of the errors in chirp mass and 
77 for different PN orders. There is significant improvement in the estimation by the use of the 3.5PN phasing instead 
of the 2PN one especially for more massive systems. For a prototypical system of a binary BH each of mass IO^Mq, 
we find that the chirp mass and rj improve by 11% and 39% respectively. Improvement is higher for more massive 
systems. For a 2 x lO^M© binary, the chirp mass and 77 improves by 14% and 45%. They are similar to the results 
for the ground based detectors as discussed in |35j but at an entirely different mass range. For a typical binary in 
its sensitivity band, LISA will be able to measure chirp mass with an incredibly small fractional accuracy of ~ 10^® 
and ry by about 10~^. 

Variation with mass: The estimation of the chirp mass worsens with increase in total mass of the binary whereas 
the estimation of 77 improves initially and then decreases. These effects can be understood as follows. When the 
total mass increases there are two competing effects in action: the increase in errors with mass, since signal lasts 
for smaller duration, and the variation of SNR with mass, which is a characteristic of the noise curve. For chirp 
mass the errors increase so rapidly that the variation in SNR does not affect the trend and the errors continue to 
increase monotonically with mass. For r/, there is trade-off between these two competing effects which accounts for 
the minimum in the curve. 



C. Errors in coalescence time 



Measuring the time of coalescence of a binary system is important to carry out electro-magnetic observation of the 
event associated with the binary merger. We discuss the trends in estimation of tc below. 

Fig. 121 displays the variation of errors in tc with increasing mass of the source and across the PN orders. The errors 
in tc show trends similar to that of |35l |. i.e., with increase in PN order the errors oscillate in a sense opposite to the 
mass parameters and in going from 2PN to 3.5PN there is a net degr adation in its estimation, which is about 43% 
for the 2 x 10^ Mq system considered. This was explained in Ref |33 based on the correlations between tc, M. and 
rj. It was noticed that both cmi^ ^-nd Cr^t^ are positive and follow the same trend as the error in tc. Increase in these 
correlations implied a worsened estimation of tc. 



D. Post- Newtonian convergence in the parameter estimation context 

Since the PN series is an asymptotic series, the rate of convergence of the results is a very important issue for 
detection as well as parameter estimation. We use the word 'convergence' to mean that the difference (in errors) 
between two consecutive PN orders is smaller as we go to higher orders considered. As remarked in Refs. |3lLl33 |. if 
the parameter estimation scheme is based on a lower order (2PN) phasing, the systematic errors due to the absence of 
higher order terms may be more than the statistical errors caused by the noise. Since we study here the implications 
of 3.5PN phasing, we examine the convergence of the series based on our results for different PN orders. 

As Fig. n reveals, though there will be improvement by using the 3.5PN phasing instead of the 2PN one, most of 
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FIG. 2: Variation of errors in tc with the total observed mass for different PN order restricted waveforms for LISA. Pattern 
averaged waveform is used. The convergence of the results is evident from the plot. Sources are assumed to be at 3 Gpc. 



TABLE I: Variation of errors in different parameters and number of GW cycles with PN order. Errors are calculated with the 
pattern averaged waveform. The system considered is a binary of mass 2 x IO^Mq at a luminosity distance of 3 Gpc. Most of 
the features regarding the improvement in parameter estimation, convergence of the PN series and correlation with number of 
GW cycles are captured in this table. 



PN Order 


Ate 


AM/M 


A77/7? 


cycles 




(sec) 


(io-«) 


(10-*) 




IPN 


0.2474 


6.217 


6.287 


2414.03 


1.5PN 


0.3149 


3.648 


1.427 


2310.26 


2PN 


0.3074 


3.694 


1.572 


2305.52 


2.5PN 


0.3947 


3.320 


0.9882 


2314.48 


3PN 


0.3435 


3.377 


1.033 


2308.73 


3.5PN 


0.4399 


3.300 


0.9661 


2308.13 



the improvement seems to come from the transition between 2PN and 2.5PN after which the series continues to show 
its characteristic oscihatory behaviour, but with smaller magnitude, suggesting that phasing at orders higher than 
3.5PN may not cause much improvement (see Table I). 



E. Parameter estimation and Number of GW cycles 

In Ref . [S^ , the correlation between the improvement in errors across different PN orders and the number of total 
and useful GW cycles [s^ was studied. It was found that though they are good indicators of how the errors at each 
order vary with the total mass of the system, they alone cannot explain the variation of errors across different PN 
orders in the context of ground based detectors. We confirm this feature in the LISA context. We recover the results 
in Table I and II of Ref as a check of our calculation. Table I shows how the errors and total GW cycles vary 
with increasing PN orders. One would expect an improved (worsened) estimation if number of GW cycles increases 
(decreases) between two consecutive PN orders if they were solely responsible for the trends. From the table it is clear 
that from IPN to 1.5PN and 1.5PN to 2PN, the errors in chirp mass and 77 do not conform to the above expectation. 
The same is the case in going from 3PN to 3.5PN. Finally, trends in tc being opposite to that of the other two mass 
parameters lead one to conclude once again that the number of GW cycles is not sufficient to understand the variation 
of errors with PN orders. 

Irrespective of whether the total number of GW cycles is very high (~ 10^) (as in the case of LISA) or low (~ 
hundreds) (as for the ground based detectors) the PN trends in parameter estimation arc too complicated to be 
explained solely in terms of this. However the variation of errors with mass can be understood on the basis of number 
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of the cycles and variation in SNR (again, not in terms of only one of them). 

F. Parameter estimation for unequal mass binaries 

Lastly we perform a similar analysis for unequal mass systems where rj < 0.25. For larger mass ratios (77 ~ 10^^ or 
smaller in our case), the Fisher matrix becomes ill-conditioned j6|| and hence wc restrict ourselves to the inspiral of 
a binary consisting of an intermediate mass BH (IMBH) and a SMBH rather than stellar mass-SMBH inspirals. An 
IMBH of 10'' Mq inspiralling into a SMBH of IO^Mq constitutes our prototypical system. This system at a distance 
of 3 Gpc will have a SNR of a few hundreds. We find that the improvement due to the inclusion of higher order terms 
is more dominant here than for the equal mass binary case. For the prototypical system considered above, wc find 
an improvement of 11% for the chirp mass and 52% for r] in the pattern-averaged case. For a 10"* — lO^M© binary, 
where SNR is ~ 100, the improvement is even more: 20% and 62% respectively^. 

This larger improvement for the unequal mass binaries is not a special feature of the LISA noise curve; for the 
ground based detectors also a similar feature exists. But unlike in the LISA case, where many such unequal mass 
binaries arc astrophysically plausible, for the ground based detectors such sources are not prototypical. 

G. Effect of low frequency cut-off chosen 

All the calculations so far, and hence the results, have been based on the optimistic possibility that the low frequency 
sensitivity of LISA can be extrapolated from lO^'' Hz to 10^^ Hz. As argued in Ref this significantly improves 
the parameter estimation. We quantify the effect of this choice of lower cut-off by comparing our previous results 
with the one where LISA is assumed to be 'blind' below 10~^Hz. 

As discussed in Sec. Ill Al we have chosen the lower limit of integration in all the calculations assuming that the 
system is observed for one year before coalescence, when / ~ f\so- By this procedure, the lower cut-off for a 2 x lO'^'^M© 
binary is about 10""* Hz. This means for systems with masses higher than 2 x (and hence a lower /iso) the 

lower limit of integration for one year observation time will be less than lO""* Hz. If we assume LISA is not sensitive to 
signals below lO"** Hz, these systems will be observed effectively for less than a year leading to a significant decrease 
in the number of GW cycles and consequent degradation in parameter estimation (for a 2 x IO^A/q binary, a lower 
cut-off of 10^'"' Hz will give 2308 GW cycles as opposed to 608 if the cut-off was lO^"' Hz). Thus for binaries whose 
masses are higher than 2 x 10'^'^Mq, the choice of lower cut-off frequency will affect the results reported here. Fig. |2| 
displays the variation of 2PN and 3.5PN errors in chirp mass and rj corresponding to the two different lower frequency 
cut-offs we have chosen. Indeed, as is evident from the plot, the errors start to deviate for binaries whose masses 
are greater than 2 x 10^-^ Mq. For a 2 x lO^M© system, using 10~^ Hz as cut-off instead of 10~* Hz improves the 
estimation of chirp mass by about 150 times and that of 77 by 40 times. These results confirm the need to push to the 
extent possible the lower frequency sensitivity of LISA. 

Regarding the improvement in parameter estimation in going from 2PN to 3.5PN, calculation with a cut-off of 10~* 
Hz shows that for a 2 x IO^Mq binary the difference in going from 2PN to 3.5PN would be 22% and 60% (as opposed 
to 13% and 45% with 10~^ Hz) for chirp mass and 77. The number of GW cycles for a cut-off of 10~^ Hz is just 7 
whereas with a 10~^ Hz cut-off it is about 540. Therefore with a cut-off of 10~^ Hz, one is only observing the very 
late inspiral of the system'' whereas with a cut-off of 10~^ Hz, the inspiral phase is dominant. The significantly larger 
variation in parameter estimation in going from 2PN to 3.5PN with the 10"'' Hz cut-off could be due to the generally 
accepted fact that higher order terms in the phasing formula are more important as one approaches the last stable 
orbit. 



IV. PARAMETER ESTIMATION WITHOUT PATTERN AVERAGING 

Having discussed in detail the various aspects of parameter estimation using pattern averaged waveform and factors 
which affect the process, we now turn our attention to the parameter estimation without pattern averaging. Obviously 



^ Using a calculation of the number of GW cycles, Ref l57l has emphasized the need for higher order PN modelling of the IMBH-SMBH 
binaries. The effects of eccentricity could also play an important role in the dynamics of such binaries. 

* Inspiral waveforms would be inadequate in this studies. Theoretical approximants, e.g. Effective one body l58l l59l may have to be 
employed to model this phase of the binary's dynamics. 
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FIG. 3: Variation of errors of chirp mass and 77 with the total observed mass for different PN order restricted waveforms for 
LISA. Pattern averaged waveform is used, for the two choices of the lower frequency cut-offs /i — 10^'' Hz and /2 = 10^^ Hz. 
Sources are assumed to be at 3 Gpc. 



the source's location, orientation and luminosity distance to the source gets added to the space of parameters which was 
four dimensional in the previous case. Unlike the ground based detectors, LISA can measure the distance, location 
and orientation of the source with a single detector because of the modulations due to its orbital motion l37j . 
Besides, using LISA as a two detector network improves the estimation of angular resolution of the source |33|. In this 
section we will discuss the improvement brought in by the higher order terms using non-pattern averaged waveform 
for LISA. We check our code, which now includes the pattern functions, by reproducing the results of [33 at 1.5PN 
with their signal and noise models. 

But as mentioned earlier the strong dependence, as we shall discuss, of the errors on the angular variables makes 
our analysis, which is for selected values of them following Cutler [3 71. le ss rigorous. The best way to deal with 
this situation is to perform Monte Carlo simulations, similar to [1,SUE3- However we notice that certain general 
conclusions can still be drawn from our limited analysis which is the topic of discussion of this section. 



A. Comparison of different detector configurations 

We compare the estimation of errors with different detector configurations now. The final errors, in comparison 
with the pattern averaged case, will depend on (i) the value of SNR corresponding to the set of angles chosen (ii) The 
change in SNR relative to the pattern averaged case and (iii) worsening of errors due to the introduction of the new 
parameters. 

For the pattern averaged case, since there are no pattern dependent parameters, the total parameter space is 
essentially 4 dimensional, the parameters being {^c, 0c, A^, '7}- When we make a transition from the pattern averaged 
to the non-pattern averaged case, five new parameters {Dl, Mi, Ms, i^i, 1^5} corresponding to the distance to the 
source, its location and orientation are added to the space of parameters. The parameter space is now nine dimensional 
significantly higher than the earlier four dimensional one. This increase in dimensionality of the parameter space leads 
to an increase in the errors of the four existing parameters. On the other hand, the introduction of pattern functions 
results in a change in SNR depending on the four angles chosen. The final picture is a complex interplay of all these 
features. 

The errors in chirp mass and 77 together with SNR is displayed in Fig. ^ for three set of angles from the seven given 
in [3^ . The three configurations corresponding to the pattern averaged, non-pattern averaged with one detector 
and the two detector network (without pattern averaging) are considered. As is evident from the plot, the SNR and 
hence the errors crucially depend on the location and orientation of the source. Also there can be orientations which 
may have lower SNR than for the pattern averaged case (see third column e.g.). However an interesting point from 
the figure (and from those runs corresponding to other values of angles which are not displayed) is that among the 
three effects which influence the parameter estimation without pattern averaging, the value of SNR seems to be the 
dominant one. Smallest errors in the plot correspond to the configuration with the largest SNR. Between the pattern 
averaged and the one detector case, the effect of correlations due to the additional parameters play a significant role. 
For example, in the second column, though one detector case has larger SNR, the errors are smaller for the pattern 
averaged case for this reason. We do not discuss the percentage improvement arising from the higher order PN terms 
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FIG. 4: Variation of the signal to noise ratio, relative errors in chirp mass and relative errors in r\ with the total observed mass 
of the binary for different choices of location and orientation of the source. Sources are assumed to be at a luminosity distance 
of 3 Gpc and a non-pattern averaged waveform is used. Angle:l corresponds to {p,L = 0.5, p,s ~ —0.8, (j)L = 3, (^s = 1}. Angie:2 
is {p,L = 0.2, fis ~ —0.6, 4>L = S,4>s = 1} and Angle:3 {fiL ~ 0.8, fls = 0.3, (j>L — 2,(f>s = 5}. The errors thus depend very much 
on the position and orientation of the source in the sky. 



since it depends very mueh on the orientation of the source. Exhaustive Monte Carlo simulations have to be performed 
to have a detailed understanding of this. 

We conclude with a remark about the estimation of angular resolution and distance as observed from the limited 
set of angles we have considered. The estimation of distance and angular resolution is not improved much because of 
the additional phasing terms. This is not surprising, since the additional terms in the phasing formula do not carry 
any information about location or orientation of the binary. One may need to go beyond the restricted waveform 
model of the waveform in order to achieve this. Some preliminary studies in this regard 1^ are 

consistent with the same. Going beyond the restricted waveform approximation would mean includin g th e amplitude 
corrections to the waveform from the two GW polarizations, currently completed up to 2.5PN order |23ll2^. This is 
because the amplitude terms are functions also of the angular positions of the source in the sky, the introduction of 
which could break different degeneracies, allowing better parameter estimation |40l l4lj . 



V. SUMMARY AND FUTURE WORKS 



The significance of higher order phasing terms is investigated in the LISA case for different sources using a pattern 
averaged waveform model. Using the 3.5PN inspiral waveform instead of the 2PN one which is currently employed in 
the GW experiments, mass parameters can be estimated with improved precision for LISA. For an equal mass binary 
of 2 X IO^Mq at a luminosity distance of 3 Gpc, the improvement in chirp mass is ~ 11% and that of is '--^ 39%. 
The FN series shows convergent behaviour beyond 2.5PN order. The number of GW cycles is a good indicator of 
how the errors vary with mass but not across different FN orders. The improvement in parameter estimation is more 
pronounced for binaries with unequal masses. The problem of unequal mass binaries should be revisited with the 
inclusion of spin effects to higher orders as spin effects play a dominant role in the dynamics of such systems. Inclusion 
of the spin effects at 2.5FN order in the phase [6^l63 | would help us go beyond the results of 0,113 for the spinning 
case and will be exciting especially for the equal mass case. The effect of the lower cut-off frequency we have chosen 
(10~^ Hz) on the parameter estimation is studied by comparing the calculation with the more modest cut-off of 10~^ 
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Hz. The estimation of source location and orientation is also studied using the non-pattern averaged waveform for 
selected source directions and orientations. They do not improve significantly by the use of the restricted 3.5PN 
template. One may need to go beyond the restricted waveform approximation in order to achieve it. Especially, 
implications of the full waveform with the 2.5PN polarizations of Ref together with the 3.5PN phasing of 
would be an interesting exercise to carry out. 

Exhaustive Monte Carlo simulations are required in order to understand the effect of higher FN terms for LISA 
when the non-pattern averaged waveform is used since the parameter estimation in this case strongly depends on the 
location and orientation of the source. We plan to take up this problem in the near future. 

It may be interesting to employ Markov chain Monte Carlo methods for parameter estimation [l^, and compare 
the results with the ones obtained using the covariance matrix. Extending this to the present context would surely be 
interesting and should be addressed. Finally, the parameter estimation for the SMBH inspirals may eventually have 
to be done using the time delay interferometry variables (see [63 | and references therein). 

The very high accuracy parameter extraction possible with LISA will make it a useful tool of astrophysics and 
provide thorough probes of strong field aspects of gravity in the future. 
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